Spin Nernst effect in the absence of a magnetic field 
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We study the spin Nernst effect of a mesoscopic four-terminal cross-bar device with the Rashba 
spin-orbit interaction (SOI) in the absence of a magnetic field. The interplay between the spin Nernst 
effect and the seebeck coefficient is investigated for a wide range of the Rashba SOI. When no peaks 
appeared in the seebeck coefficient, an oscillatory spin Nernst effect still occurs. In addition, the 
disorder effect on the spin Nernst effect is also studied. We find that the spin Nernst effect can be 
enhanced up to three-fold by disorder. Besides, due to the interface effect, the counter-propagating 
of the charge current to the direction of the temperature gradient is possible for a nonuniform 
system. 
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With the development of the micro-fabrication technol- 
ogy and the low-temperature measurement technology, a 
great amount of efforts have been paid for the research of 
the thermoelectric properties in the last two decades^. 
Comparing to the conductance, the thermoelectric co- 
efficients of electronic systems are more sensitive to the 
details of the density of states^—, which is very important 
for the design of the electronic devices. The thermopower 
(seebeck coefficient) of the quantum dot was measured in 
the last few years^. Recently, the Nernst effect, a Hall- 
like thermal effect, has been theoretically studied 7 and 
had been detected, for example, in bismuth 8 in which, 
with the existence of a perpendicular magnetic field, a 
transverse current is induced by the longitudinal ther- 
mal gradient. 

In the spintronics area, the spin thermal coefficients 
are also of focus recently^—. In a recent paper, by con- 
sidering a system with a spin-orbit interaction (SOI), the 
Nernst effect and a novel thermal effect, the spin Nernst 
effect, have been fully studied in a two-dimensional elec- 
tron gas^. It is found that, because of a perpendicular 
magnetic field B, the Nernst signal exhibits a series of 
peaks. When the SOI exists, the peaks split and the spin 
Nernst effect appears. With a small B or a large SOI, 
the spin Nernst effect becomes more pronounced. It also 
shows that the spin Nernst effect is easier to be affected 
by disorder than the Nernst effect. 

There is no doubt that a perpendicular magnetic field 
B is essential for the existence of the Nernst effect. How- 
ever, in the spin Hall effect, the transverse spin current is 
due to a SOI rather than a perpendicular magnetic field. 
Similarly, for the spin Nernst effect, B may not be needed 
either. One may suspect that the spin Nernst effect is in 
fact the combination of the existence of thermopower and 
a SOI. Thus, the focus of the current work is to study 
the spin Nernst effect in the absence of a perpendicular 
magnetic field, and its interplay with the thermopower. 

In this paper, the property of spin Nernst effect is de- 
veloped in a two-dimensional electron gas system with a 
Rashba SOI but without a perpendicular magnetic field 
B. For this set-up, the Nernst effect disappears thus we 



focus on the spin Nernst effect - a transverse spin cur- 
rent induced by a longitudinal thermal gradient A7~. A 
traditional way to analyze such a Hall-like system is to 
add vertical probes to detect the transverse properties. 
Thus we set a four-terminal cross-bar sample, as shown 
in Fig.l— . A longitudinal thermal gradient AT is added 
between the leads 1 and 3. This thermal gradient induces 
a transverse spin current J s in the closed boundary condi- 
tion with a SOI, which can be measured at leads 2 and 4. 
The seebeck coefficient of such a system can be directly 
measured at leads 1 and 3. 



By using a tight-binding model and the Landauer- 
Buttiker (LB) formula, the spin Nernst coefficient N s 
(N s = J S /AT) and the seebeck coefficient S (S = 
—Ay/ AT) are calculated. The Rashba SOI used in our 
calculations covers a wide range with some beyond the 
accessibility of today's sample. The seebeck coefficient S 
shows a few peaks consequently when the fermi energy 
Ep goes through the energy band. Due to the interface 
of our setting (zero Rashba SOI at lead 2,4), we find a 
negative S. It is confirmed that spin Nernst effect can 
not be simply thought as the combination of the seebeck 
coefficient and the Spin hall effect^. A big spin Nernst 
coefficient N s can be found with a zero seebeck coeffi- 
cient S. However, when the peaks of seebeck coefficient 
occur with a non-zero Rashba SOI, the spin Nernst ef- 
fect exhibits big amplitude or sometimes also peaks. The 
Fermi energy Ep also affects N s . When the Fermi energy 
Ep is close to the bottom of the energy band (— At), the 
oscillatory amplitude of N s becomes more pronounced. 
The effect of disorder on N s is also investigated. When 
Ep = — 3.8t, we can see a large increase of N s with in- 
creasing of the strength of disorder. Its value at the peak 
is about three- fold of that without disorder. In addition, 
we find that the strength of disorder when N s vanishes, 
indicating that the system goes into an insulating regime, 
is independent of the Fermi energy. 



In the tight-binding representation, the Hamiltonian 
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FIG. 1: Schematic diagram of the four-terminal cross-bar 
sample. The area with SOI is marked by gray. A thermal 
gradient AT is applied between the longitudinal lead-1 and 
lead-3. 



with SOI can be written as>^, 

+4+s x .A- tl + i°yVn)vv'C ia , + H.c] (1) 



where cJ _(ci CT ) is the creation (annihilation) operator of 
electrons in the site i = (n, m) with spin a, and 5x and Sy 
are the unit vectors along the x and y directions. e\ is the 
on-site energy, which is set to everywhere for the clean 
system. When the center region is a disorder system, £i is 
set by a uniform random distribution [-W/2,W/2]. Here 
t = h 2 /(2m*a 2 ) is the hopping matrix element with the 
lattice constant a, I is a two-dimensional identity matrix. 
The strength of Rashba SOI is represents by Vr = ah/ 2a, 
where a is the Rashba spin-orbital coupling. Vr is set to 
zero in the lead-2 and lead-4. 

Considering a small temperature gradient AT on the 
longitudinal lead-1, 3, we can set the temperatures 7i — 
T + AT/2, %=T- AT/2, T 2 = T = T. The charge 
current in lcad-p can be written as J pe = e(I p -\-\- I p {) and 
the spin current is J ps — (H/2)(I p f — I P i). Here I pa is 
the particle current in lead-p with a equals to t or J.. I pa 
can be obtained by the LB formula>2iiS 
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dE T pa , q (E){f p (E) - f q (E)} (2) 



where T pa>q (E) is the transmission coefficient from the 
lead-q to the lead-p with spin a and E is the energy 
of the incident electron. f P (E) is the electronic Fermi 
distribution function of the lead-p. 

The spin Hall current in lead-2 and lead-4 can be calcu- 
lated with the closed boundary condition in both lead-1, 3 
and lead-2,4, i.e. V x = V 3 = and V 2 = V 4 = 0. From 
symmetry of the system, we know that J 2s = —Ji S 9 - 
After the Taylor expansion, the spin Nernst coefficient 
N s = J 2s /AT can be reduced to: 



N s = 



= ^J dE(AT 23 - AT 2 i) 



E — Ep 
k B T 2 




FIG. 2: N s (red solid) and S (black dotted) vs. Fermi energy 
Ef for different Rashba Vr. The (scaled) transmission coef- 
ficient 21,2 + Tl,3 (thin blue dashed) and spin transmission 
coefficient AT2,3 (thin blue solid) are also shown. The other 
parameters are T = 0.01, and L = 19a. 



here AT 2p — T 2 ^^ p — T 2 ^ p , and / is the zero order of Tay- 
lor expansion of the Fermi distribution function, it is the 
same for all four leads, f(E) — I / {exp[(E — Ef) / ksT] + 
!}• 

For the calculation of the longitudinal seebeck coeffi- 
cient S, we need the open boundary condition at lead-1, 3, 
i.e. Ji e = Jze = to find the difference AV^ = V\ — V3. 
Different from a quasi-one-dimensional 2-leads system^, 
the extra leads-2,4 also affects the longitudinal seebeck 
coefficient S of the entire system. For example, with 
a perpendicular magnetic field B, the longitudinal see- 
beck coefficient S is affected by the bias in leads-2,4, 
V 2 and V4. However, without B, the sample's symme- 
try increases from C 2 symmetry to D 2 symmetry, i.e., 
we have T 1)2 = T M , Here T lj2 = 7if,2+ ?u,2- After the 
Taylor expansion, we can get the longitudinal seebeck 
coefficient S = -AV/AT as: 



S 



1 / dE (T x , 2 + T 1A + 2Ti, 3 ) (E - E F )f (1 - /) 



T 



JdE(T 



1,2 + T 1A + 2Ti, 3 ) / (1 - /) 



/(I-/), (3) 



( 4 ) 

The equation above shows that, even with a higher sym- 
metry, the longitudinal seebeck coefficient S is still af- 
fected by the transport properties from lead-2 and lead- 
4. 
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FIG. 3: (a) A simple model: Current J because of voltage 
gradient (red dashed) and thermal power S because of ther- 
mal gradient (blue dotted) vs. Fermi energy Ef at a two-lead 
system with Rashba Vr = 0. The (scaled) transmission func- 
tion Ti t 3 is also shown (black solid). The plot in the small box 
shows ft — fa with temperature difference, (b) and (c): the 
eigen energy of the lead E n> k v.s. longitudinal wave vector 
k y (units: l/2a) for different Vr 



With the £>2 symmetry, the relationship between S 
and N s can be further derived. In fact, we can rewritten 
S = (A^ + A?) I [J dEF (e) (at + a 1 )] . The D 2 sym- 
metry gives AT 23 = -AT 2 i. Noticing T 3t4 = T 3 ^i, 
the spin Nernst coefficient can be simplified as N s — 
(A± - At) / (2irk B T). Here e = E - E F and F (e) = 
/ (1 — /), at denotes the spin up term: at = TV^i+Ta-^i, 
and at the spin down term at = T^i + T^i, we also use 
the notation of the integral term At — j dEsF (e) at j T 
and A^ = J dEeF (e) at / T. Because of the symmetry, 
only leads- 1,2, 3 are used in the simplified expression of 
S and N 8 , we only need the upper half of the sample for 
our investigation. In fact, at (a^) and A^ (A^) reflects 
transport properties of spin-up (spin-down) electrons in 
the upper half of the sample. Roughly speaking, S can 
be seen as the sum of spin-up and spin-down terms, while 
N s as the difference of them. 

In the numerical calculations, t = fi 2 /(2m*a 2 ) is set 
as the energy unit. If taking the effective electron mass 
to* = 0.05TO e and the lattice constant a = 12.5nm, t 
is about bmeV . Temperature is fixed by k B T = O.Oli, 
which is about IK. The size of center region is L = 19a, 
about 237rtTO. In a reasonable experimental range thus 
far Vr £ [0, 0.l}^ 2 -. However, in order to thoroughly study 
the relationship between the spin Nernst coefficient N s 
and the seebeck coefficient S, we extend the range of 
Vr up to [0, 1] in our calculation. 

Fig. 2 shows the spin Nernst coefficient N s and the see- 
beck coefficient S versus the Fermi Energy Ep in the 
clean system (W — 0). It is clearly seen that the see- 



FIG. 4: (N s (red solid) and S (black dashed) vs. Rashba SOI 
Vr for fermi energy Ef = —3. For compare, the blue dotted 
line shows the seebeck when lead-2,4 have the same Vr as 
lead-1,3. The other parameters are T = 0.01, and L = 19a. 



beck coefficient S peaks at the positions where there 
are step-changes of transmission function T\p, + T\^. 
These peaks can be explained by a simple model only 
with a 2-lead system without the Rashba SOI, shown in 
Fig. 3(a). The transmission coefficient is a step func- 
tion (solid-black curve). The reason is as follows. The 
sample can be considered as a multi-channel system at 
a low temperature (here T ~ IK). When fermi energy 
increases, more channels in the lead are used to trans- 
port current. Thus, AV of two leads as well as the cur- 
rent increase with increasing of fermi energy (red-dashed 
curve). However, the S (blue-dotted curve) can not ac- 
cumulates while Ep increases, it only peaks while the 
channel number changes and S is close to zero with a 
fixed channel number. This can be seen from the LB 
formula l|2~|). if lead-p and lead-q have different tempera- 
tures, f p (E, T+AT)—f q (E,T — AT) is an antisymmetry 
function of E — Ep (see plot in small box of Fig. 3(a)): 
when E < Ep, f p < f q , current flows from lower tem- 
perature lead to higher temperature one; when E > Ep, 
current flows in the opposite direction. Only when the 
two flows are not equal, i.e. T V a,q has an antisymmetry 
part, we can have a nonzero current. Thus for Fig. 3 (a), 
only when T13 is at the step-change point, it has anti- 
symmetry part and can give a non-zero S. 

This conclusion can also be used to analyst spin- 
involved quantities. From Fig. 2, we can see that N s (red 
solid line) shows an oscillatory structure. Besides the 
peaks at Vr = (N s is zero at this point), the mag- 
nitude of N s oscillation is also large at the peaks of S; 
but at the exact maximum point of S, where Vr is quite 
small (Vr < 0.1), N s is generally close to zero. This is 
because the spin transmission coefficient AT2.3 generally 
has an extreme value when the transmission coefficient 
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FIG. 5: (a). 7V S vs. the strength of disorder W for different 
Fermi level Ef = — 3.8(solid blue), Ef = — 3.6(dashed red), 
E F = -3.0(dotted black), E F = -2.2(dotted dashed green); 
(b). N s vs. Fermi energy Ef = —3.8 for W — (dashed 
black) and W = 1.4 (solid black). Other parameters are Vr = 
0.05, T = 0.01, and L = 19a. 



jumps at a step. Around an extreme value, any func- 
tion is almost symmetry, thus one only can get a low 
value of N s . While at both sides of the extreme value, 
AT2 i 3 monotonically increases or decreases, we can get 
a local maximum magnitude of N s . Now why AT 2 ,3 has 
an extreme value at a peak of S for a small Vr. Due 
to the Rashba SOI, each eigen-energy band splits into 
two sub-bands with opposite spin directions. These two 
sub-bands degenerate at k y = 0, and the lower sub-band 
has two valleys below this degenerate point. The two 
sub-bands are very close to each other when Vr is small. 
If the lower sub-band of high level (for example, E\ t k y ) 
has the similar spin direction with the upper sub-band of 
low level energy (for example, Eq^ ), AT2 ; 3 continually 
increases when Ep goes from the upper band of £ufc y 
to the two valleys of lower sub-band of E\^ , and than 
rapidly decreases when Ep goes through the degenerate 
point {k y = 0) of Ei.k y , thus we get a peak in AT2.3; 
otherwise we get a valley in AT2.3- 

When Vr is very big, we can see the external peaks 
for both N s and S. For example Vr = 0.45 in Fig.2(d), 
close to the 2 nd and 3 rd main peaks of S, we can see 
a very sharp sub-peak of N s . In fact, these are also the 
small peaks of S, though not very big. This is because for 
these two band (see Fig. 3(c)), the two valleys of the lower 
sub-band is far from the degenerate point at k y — 0, the 
two channel of these two sub-bands is separated. Thus 
we can see two peaks. At this time, the change of spin 
transmission coefficients can be roughly thought as the 
change of transmission coefficients, thus we can see N s 
peaks at the S's peak. 

In Fig. 4, we show the spin Nernst coefficient N s and 



the seebeck coefficient S versus the Rashba SOI Vr in 
the clean system (W — 0) for Ep = —3. The seebeck 
coefficient S decreases and maintains for a small value for 
quite a while before shows another peak. This is because 
increasing Vr moves the energy bands and makes them 
go through the fermi energy. It should be mentioned 
that we found the negative seebeck coefficient S (Fig. 
2d), which means a longitudinal current occurs in the 
opposite direction of the temperature gradient AT. This 
is due to the boundary conditions Vr = at leads-2,4. 
As a compare, we also show S for a uniform system, 
i.e. leads-2,4 having the same strength of Vr as in the 
sample. For this situation, the seebeck coefficient S is no 
longer negative. In fact, when the Rashba SOI is absent 
in the leads-2,4, an interface between Vr — and Vr ^ 
ocurrs^, this interface causes additional scattering for an 
incident electron. In some special case like Ep = —3, 
this may make the electrons below Ep easier to transport 
than the electrons above Ep, thus a negative S. 

Finally we discuss the disorder effect on the spin Nernst 
effect. Fig. 5 shows N s versus disorder strength W for dif- 
ferent Fermi energies. The calculations are averaged over 
500 disorder configurations. Around W < 1.7, N s shows 
an oscillatory structure. N s changes sign with increas- 
ing of the disorder strength (see Ep = —3.6 and —3.0 in 
Fig. 5). It is interesting to see that, comparing to a clean 
system (W = 0), N s can be unexpectedly increased by 
disorder W. This is because the disorder changes the os- 
cillating structure of N s (see Fig. 5b). As expected, the 
disorder decreases the strength of oscillating, however, it 
also shifts the peak positions of N s . It is possible to have 
a peak in N s at finite disorder while it is almost zero 
initially at clean limit. In Fig. 5a, around W = 1 ~ 1.5, 
for the Fermi level Ep = —3.8, —2.2, we can see that N s 
is up to about three times of N s at W — 0. The behav- 
ior of N s v.s. W is very apparent when the Fermi level 
Ep is close to the bottom of energy band (Ep = —4). 
For Ep = —3.8, N s is much bigger than those at other 
Fermi levels, and we can see a very remarkable peak at 
about W = 1.4. For E F = -2.2, N s begins from -0.005, 
changes its sign at about W ~ 1.25 and than increases, 
again reaches to 0.005 at about W ~ 1.75. With a very 
big disorder, N s should go to zero as system enters into 
an insulating regime. We find that the zero of N s occurs 
at W = 3 for Vr = 0.05. This is roughly independent of 
the locations of the Fermi energy. 

In summary, in the absence of a perpendicular mag- 
netic field, the interplay between the spin Nernst ef- 
fect and the seebeck effect is investigated in a two- 
dimensional cross-bar with a spin-orbit interaction. The 
spin Nernst effect exhibits an oscillatory structure for a 
wide range of the Rashba SOI. With a large Rashba SOI, 
the N s oscillation has a peak when the seebeck coeffi- 
cient possesses one. However, the inverse condition is 
not always satisfied, namely, the seebeck coefficient can 
be almost zero while N s has a peak. The disorder ef- 
fect on the spin Nernst effects is also studied. We find 
that disorder can enhance N s up to three times for some 
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Fermi levels. In addition, the disorder can also change 
the sign of spin Nernst effect. Moreover, the limit of dis- 
order where N s goes to zero is independent of the Fermi 
energy. 
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